Magnetoelectric coupling at the interface of BiFe03/Lao.7Sro.3Mn03 multilayers 



M.J. Calderon/ S. Liang,^'^ r. Yu,"^ J. Salafranca,^ S. Dong,^ S. Yunoki/'^ 

L. Brey,^ A. Moreo,^'^ and E. Dagotto^'^ 
^ Instituto de Ciencia de Materiales de Madrid (CSIC), Cantohlanco, 28049 Madrid, Spain 
^Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, USA 
^Materials Science and Technology Division, Oak Ridge National Lab, Oak Ridge, Tennessee 37831, USA 
Department of Physics and Astronomy, Rice University, Houston, Texas 7705, USA 
^ Departamento de Fisica Aplicada III, Universidad Complutense de Madrid, 28040 Madrid, Spain 

^Department of Physics, Southeast University, Nanjing 211189, China 
'^Computational Condensed Matter Physics Laboratory, RIKEN, Wako, Saitama 351-0198, Japan 
^ CREST, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 

(Dated: December 8, 2010) 

Electric- field controlled exchange bias in a heterostructure composed of the ferromagnetic mangan- 
ite Lao.7Sro.3MO3 and the ferroelectric antiferromagnetic BiFe03 has recently been demonstrated 
experimentally. By means of a microscopic model Hamiltonian we provide a possible explanation of 
the origin of this magnetoelectric coupling. We find, in agreement with experimental results, a net 
ferromagnetic moment at the BiFe03 interface. The induced ferromagnetic moment is the result 
of the competition between the Cg-electrons double exchange and the t2g-spins antiferromagnetic 
superexchange that dominate in bulk BiFe03. The balance of these simultaneous ferromagnetic and 
antiferromagnetic tendencies is strongly affected by the interfacial electronic charge density which, 
in turn, can be controlled by the BiFe03 ferroelectric polarization. 

PACS numbers: 75.47.Gk, 75.10.-b, 75.30.Kz, 75.50.Ee 



Introduction. The quest for efficient electric field con- 
trol of magnetic properties has encouraged research on 
materials with a strong coupling between the magnetic 
and dielectric degrees of freedom JT^. Such a control 
would find applications in magnetic field storage and 
sensors, constituting a major step forward in the field 
of spintronics. However, so far no bulk material seems 
to possess the required characteristics, including working 
near room temperature. For these reasons, the field ef- 
fect device presented in Ref. [7 signals a new route [3 to 
achieving those goals by growing a few nanometers thick 
layer of Lao.rSro.sMnOa (LSMO), a ferromagnetic (FM) 
metal, on an antiferromagnetic (AFM) and ferroelectric 
(EE) material, BiEeOa (BEO) 0^. These experiments 
provided evidence for an induced EM moment in BEO at 
the interface. This magnetic moment is strongly affected 
by the BEO polarization, which results in an electric field 
control of the LSMO exchange bias (EB), and concomi- 
tant control of the LSMO magnetization. 

Magnetic moments induced at the interface of 
perovskite-based oxide materials have been previously 
reported. Eor example, a net magnetic moment was in- 
duced in an AEM manganite when grown in a multi- 
layer with a EM manganite [5l [6] , at the interface of a 
superconducting cuprate with Lao.67Cao.33M03 [7 , and 
at LSMO/SrTiOa interfaces |8j. Electronic energy loss 
measurements indicate that a charge redistribution takes 
place [71 E] . Theoretically, the origin of the induced mag- 
netic moment in [5^ was explained in terms of charge 
transfer and a double exchange (DE) type interaction [9 . 
Charge transfer, together with orbital reconstruction, is 
also believed to play a role at cuprate/Lao.67Cao.33M03 
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FIG. 1. (color online) Proposed spin order in the BFO-LSMO 
heterostructure (only half of the structure is shown, the other 
half is symmetric). The two first (1,2) and two last (11, 12) 
planes have ionic charge Z=+2, and correspond to BEG. Lay- 
ers 3 and 10 define the interface, with Z=(0.7 + 2)/2. The 
central planes (4 to 9) have Z=0.7 and correspond to LSMG. 

interfaces [7] . In [H [10] the BEO induced EM moment in 
LSMO/BEO heterostructures was attributed to Ee-Mn 
hybridization, which is associated with charge transfer. 
The recently observed coupling of the EB with an in- 
duced magnetization near the interface, simultaneously 
controlled by the EE polarization [2] |4], defines a new 
complex phenomenon that requires a better theoretical 
understanding [11 . 

In this letter, a microscopic model for the transition 
metal <i-electrons is shown to explain the main proper- 
ties of the BEO/LSMO interface. Within this model, 
the magnetic moments of the Ee ions develop a net 
EM moment close to the interface as a consequence of 
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both charge and orbital redistribution, while they remain 
AFM ordered far from the interface. Interestingly, the 
direction and magnitude of the BFO moment with re- 
spect to the magnetization of the LSMO layer depend on 
the charge density at the interface. Small changes in the 
charge density due to the switching of the FE polarization 
in BFO produces large modifications in the direction of 
the magnetic moment induced in the Fe ions. This leads 
to the experimentally observed EB [2]. Our scenario is 
qualitatively different from that proposed in [4] , which is 
based on the first LSMO layer, and that of [TT], which 
requires a spin-orbit coupling. 

Model Fig. [l] sketches the 4 x 4 x 12 superceh used in 
our calculations. Bulk LSMO (BFO) is in a FM (G-type 
AFM) state. A possible magnetic reconstruction near 
the interface is also indicated. Departures from the cubic 
perovskite lattice, as induced by strain and ferroelastic 
effects, are ignored for simplicity. The z-axis size (8 layers 
for LSMO and 4 for BFO) is selected such that bulk 
behavior is recovered at the center of the composing films. 

Magnetism in the heterostructure arises from the tran- 
sition metal d electrons. In manganites, the three t2g 
electrons are localized and are approximated by a sin- 
gle classical spin. AFM superexchange with neighboring 
sites is introduced via a Heisenberg interaction [12 . The 
quantum itinerant eg electrons are described by the DE 
model in the infinite Hund's coupling limit. Several as- 
pects of the rich physics of bulk p!2HT5] and heterostruc- 
tured manganites ^ [TBHTS] have been successfully ad- 
dressed within this approach. However, model Hamilto- 
nian approaches have not been used before for BFO. 

To properly consider the effects of charge leakage and 
orbital hybridization at the interface [10] , the eg electrons 
of both BFO and LSMO must be treated on equal footing. 
The model is: 

+ E(^i + ^i-MK- (1) 

i 

Here, represents the t2g spin at site located either 
at the LSMO or the BFO side of the heterostructure, 
while Jfj^^ is the AFM superexchange parameter, cj^ 
creates an electron on an orbital centered at the tran- 
sition metal site with eg symmetry: Of, /3 = |3z^ — ^^), 
\x'^ —y'^). The hopping term is modulated by the DE fac- 
tor ^lij [12^ that depends on the angle between and Sj 
such that it is maximum for parallel alignment of spins 
at neighboring sites and zero for antiparallel alignment. 
Hopping also depends on the overlap between the a and 
/3 orbitals along the direction i — j through the geometric 
factor Of_fj [12j. In principle, the hopping parameters 
tij should depend on the material, and would be affected 
by lattice distortions near the interface. We simplify the 
calculation by assuming a uniform hopping parameter for 
the whole heterostructure {tij=t) {t is the energy unit). 
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FIG. 2. (color online) (Main panel) Energy vs. canting angle 
for different values of Jafm- (Inset) Value of the magneti- 
zation induced at the interfacial plane of BFO vs Jafm, as 
given by the position of the minima in the main panel. 

The superexchange coupling is more sensitive to changes 
in lattice constants, hence each material will be charac- 
terized by a different Jafm- 

In general, the DE term in Eq.Q favors FM config- 
urations that optimize the kinetic energy, while the su- 
perexchange term favors AFM phases. The third term 
contains the different contributions to the site potential. 
Long range Coulomb interactions are essential to control 
charge transfer across the interface. This is included in 
the Hartree approximation by setting 
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The Coulomb interaction strength is regulated by the 
parameter a, here assumed equal to 2t 0^^. For each 
x-y plane, Z as illustrated in Fig. [l] is considered, and 
the approximation is made that the background con- 
sists of point charges that occupy the transition metal 
sites. Z=2 for BFO, while Z=0.7 for LSMO. An interfa- 
cial layer is considered with an intermediate value of Z, 
to account for possible diffusions and different chemical 
environments of the transition metal ions at the inter- 
face. Vi includes the effect of the band offset between 
BFO and LSMO V^ffset^ and the surface charge density 
due to ferroelectricity Vq. The band offset is difficult to 
estimate, though expected to be small. From BFO elec- 
tronic affinity [20 and LSMO work functions [21], we 
estimate Kffset = "^lsmo - Vbfo ^ OM. However, Kffset 
is here treated as an adjustable parameter: it is first set 
to zero and the results are then checked against moderate 
changes in its value. The chemical potential /i is chosen 
so that the overall system remains charge neutral. 

Several approximations are implicit in Eq.Q. While 
the intraorbital Hubbard interaction U is effectively infi- 
nite due to the infinite Hund's coupling implicit in DE, 
the interorbital U' might be important. However, it has 
been checked that introducing a moderate U'=2t at the 
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FIG. 3. (color online) (a) Electronic charge for canting angles 
and ^^37r/2 (the dashed line is the positive charge 
background). The charge density distribution depends only 
slightly on 0. (b) Orbital order described by the expectation 
value of {Tz)=n^2_y2-n^^2_^2. The eg orbitals are equally 
populated ((r^)=0) at the BFO and LSMO bulk layers. 

mean- field level does not significantly affect the results. 
The coupling to Jahn- Teller lattice modes is not included 
either, but it is widely accepted that they do not play an 
important role in LSMO or BFO at any of the interfa- 
cial fillings of the d bands. For each set of parameters 
and different t2g spin configurations, Eq.Q was solved 
by exact diagonalization. Periodic boundary conditions 
are used, with a 5x5x1 mesh in reciprocal space. 

Results. Our first important result is sketched in Fig.[l] 
While LSMO remains FM, and most of BFO remains in 
the G-AFM state, a magnetic reconstruction takes place 
in the last atomic layer of BFO. The perpendicular ori- 
entation of the spins in bulk BFO (layer 1) relative to 
LSMO gives lower energy than a parallel orientation. The 
magnetic state in the last BFO layer can be characterized 
by a single angle 6, defined in Fig. [ij which determines 
the magnetization of BFO close to the interface, Mbfo- 
For 6>=7r/2, the last layer of BFO is FM but antipar- 
allel to the LSMO magnetization, thus Mbfo=— 1; for 
6>=37r/2, the last layer of BFO is FM and parallel to the 
LSMO magnetization, Mbfo=1; and for ^=0, 27r, there 
is no net magnetization in BFO. Regardless of the ground 
state value of ^, the main magnetic reconstruction is con- 
fined to the BFO outmost layer, since spin canting at the 
LSMO interfacial atomic layer is small. These results 
are remarkably independent of details, such as whether 
there is an interfacial layer with intermediate background 
charge. It also holds for several values of Jafm, as long as 
they are reasonable (for the proper bulk phase diagram: 
J^IT < 0-1; 0-1 < ^AFM < 0-2). 4fm'" (between 
layers 2 and 3) is expected to be some average of Jafm^ 
and Jafm- For simplicity, we use Jafm=-^afm^''''^=^afm 
and ^AFM ~ ^* 

Fig. [2] illustrates the energy as a function of for dif- 
ferent values of Jafm- At small Jafm, the energy is 



minimized at ^=37r/2: the BFO interfacial layer is FM 
and parallel to the metallic LSMO and DE dominates. 
As Jafm increases, the minimum near 37r/2 moves to- 
ward 27r. In addition, the minimum near 7r/2 decreases 
in energy and eventually has the lowest energy. In this 
case the last layer of BFO is FM and antiparallel to the 
LSMO magnetization. This trends are summarized in 
the inset of Fig. [2] where the magnetization induced in 
the last layer of BFO is plotted as a function of Jafm- 
The magnetization rotation (from 0=37t/2 to 7r/2) oc- 
curs at JAFM~0.12t and there is some small canting for 
Omt < Jafm < O.llt. For Jafm > 0.16t, the AFM 
coupling is stronger than the DE and the canting angle 
is zero, namely, the interfacial BFO layer is AFM as in 
bulk. Several results obtained in the simplified discussion 
presented here, such as Fig. |2] and others based on the 
simple assumption that the interfacial behavior is char- 
acterized by a single angle 6>, were also confirmed numer- 
ically using 4x4x8 clusters, the Poisson equation, and a 
minimization algorithm for the t2g classical spins [22 . 

The observed charge redistribution and orbital recon- 
struction confirm the importance of the DE mechanism. 
Fig.|3]shows the charge Ui and orbital polarization (r^) = 
nj;2_y2 —risz^-r^ for two values of 6. The eg-charge profile 
is mainly determined by the background charge except at 
the interface planes where there is a charge redistribution 
whose extension is controlled by the parameter a [16]: 
holes at the last atomic plane of BFO and extra electrons 
at LSMO favor a kinetic energy gain. For Mbfo=-19 the 
hopping between LSMO and the interfacial BFO layer is 
suppressed. However, even in this case the Fe (i-orbitals 
are not completely filled, and some kinetic energy gain 
occurs, mainly within the x-y plane. This asymmetry is 
evidenced by a non-zero orbital polarization (r^), at both 
sides of the interface (layers 2,3,10,11). Since there is 
no electron-lattice coupling here, the orbital polarization 
necessarily arises from asymmetries in the orbital filling 
due to kinetic energy gain. For Mbfo=1 there is hopping 
across the interface and a significant orbital polarization 
appears only at the BFO side (layers 2,11). The positive 
orbital polarization is caused by the suppression of the 
hopping in the z-direction due to the filled bands of bulk 
BFO leading to an enhancement of the hopping in the x- 
y plane. A larger kinetic energy in this case compensates 
for the Coulomb energy cost of an eg density that further 
deviates from the background value. Although the differ- 
ence in charge profile for the two angles is small, it plays 
an important role in the phenomena discussed next. 

Consider now the effect of the BFO FE polarization. 
In our study, we assume that the main effect of switch- 
ing the FE polarization is to modify the induced charge 
density at the FE surface. As a consequence, the FE na- 
ture of BFO makes the heterostructure work as a field- 
effect device [2 . The direction of the FE polarization is 
modeled by introducing an additional potential Vq at the 
surface of BFO (layers 2 and 11) which attracts (Vq < 0) 
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FIG. 4. (color online) Magnetization induced in BFO (Mbfo) as a function of Jafm and Vo for (a)yoffset = and (b)yoffs. 
(c) Electronic charge in the heterostructure for Vb = — 0.5t and Vb = 0.5t. Inset: charge at layer 2 for — It < Vb < It. 



0.6t. 



or repels (Vb > 0) the charge [see Fig. |4] (c)]. Fig. |4j 
our main result, shows the value of the induced magne- 
tization in BFO, Mbfo, as a function of Jafm and Vq^ 
for (a) zero band-offset, and (b) Kffset=0-6t. For small 
values of Jafm, a FM BFO layer appears at the inter- 
face, parallel to the magnetization on LSMO (^=37r/2). 
For larger values of Jafm and an attractive Vb, the AFM 
solution (^=0) is obtained. This is due to the fact that 
increasing the density of charge towards 2 produces a 
decrease in kinetic energy, so the gain in superexchange 
energy dominates. For large values of Jafm and repul- 
sive Vo, a FM BFO layer is obtained which is antiparallel 
to the magnetization on LSMO (6>=7r/2). This is due 
to the decrease of the charge density (away from 2) at 
the interface which produces an increase of the x-y plane 
kinetic energy, while the superexchange term gains en- 
ergy by making the spins of layers 2 and 3 antiparallel. 
Equivalent results are found for different values of the 
band-offset. Therefore, Fig. |4] explains the experimen- 
tally demonstrated control of the EB. The LSMO EB 
is determined by the magnetic order of the last layer of 
BFO, via a partial pinning to the AF BFO bulk [23], but 
this order is strongly affected by an electric field (through 
changes in Vb). A magnetic field can also influence the 
EB, by reducing the effective value of Jafm and favoring 
FM order in the last layer of BFO. 

Conclusions. A microscopic model that explains the 
recently unveiled properties of the BFO/LSMO interface 
is proposed. The charges and spins couple via the DE and 
superexchange mechanisms, and our calculations show 
that the induced magnetic moment in BFO arises from 
charge transfer between the two materials. The spin ar- 
rangement generated at the BFO interfacial layer arises 
from the frustrating effect caused by the two competing 
(FM and AFM) tendencies in adjacent layers with dif- 
ferent electronic densities [18]. Our main result is that 
changing the sign of the BFO ferroelectric polarization 
modifies the extra charge near the interface, which in 
turn strongly affects the magnitude and direction of the 
magnetic moment. This gives rise to the experimentally 



observed magnetoelectric coupling, and clarifies the ori- 
gin of the recently observed electric field controlled ex- 
change bias in LSMO/BFO heterostructures. 
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